Differential geometric and topological methods with 
MHD and plasma physics constraints 



L.C. Garcia de Andrade 1 



Abstract 



Non-solitonic examples of the application of geometrical and topological 
methods in plasma physics and magnetohydrodynamics (MHD) are given. 
The first example considers the generalization of magnetic helicity to grav- 
itational torsion loop. The second example considers also the application 
of this same torsion loop metric to some problems of helical fields in MHD 
dynamo theory. In the last example a Riemannian magnetic metric is given 
where the magnetic field itself is present in the diagonal time-independent 
metric. In this example the MHD equations are shown to be compatible to 
the geometrical Bianchi identity by making use of Cartan's differential cal- 
culus formalism. The Riemann curvature of the helical flow is also obtained. 
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1 Introduction 



In recent years several examples [1, 2, 3, 4] of the use of differential geomet- 
rical and topological methods have been presented in the literature, ranging 
from the Schief's generation of toroidal flux surfaces in MHD via soliton 
theory [1] and the hidden integrability in ideal MHD using the Pohlmeyer- 
Lund-Regge to the solar physics magnetic topology applications considered 
by Moffatt and Ricca [3] and Field and Berger [4]. In all these example the 
differential geometry of curves and surfaces [5] were used. Earlier Sivaram 
and Garcia de Andrade [6] have used non-Riemannian geometry with torsion 
to investigate the Debye problem in plasma physics. Those previous work 
using soliton theory in MHD of course made use of the Riemannian geometry 
of surfaces where the only torsion that appears was the Serret-Frenet scalar 
torsion and not the higher-dimensional [7] lie Cartan [9] torsion with the ex- 
ception of reference [6] . In this paper we strike back to the non-Riemannian 
geometry endowed with Cartan vectorial torsion to investigate torsion loops 
, previously investigated by Letelier [9] in the context of gravitational physics 
such as Einstein-Cartan gravity now also in the context of MHD, by substi- 
tuting the vector field in this metric by the magnetic field itself. We also show 
that it is possible to generalised the magnetic-helicity topological aspects of 
plasma physics can be carried over to gravitational physics. Another inter- 
esting example is the diagonal Riemannian magnetic metric. This spacetime 
metric instead of satisfying the Einstein equations ,or either Einstein-Cartan 
equations in the case of non-Riemannian loops, it fulfills the MHD equations, 
that is the reason why be called magnetic metric. The magnetic metric also 
constrain the Bianchi identity via MHD dynamo equation [10, 11]. The ideas 
applied here in the context of MHD and plasma physics keep resemblance 
with the analog gravity models which compare fluid mechanics and condensed 
matter systems such as BEC to the Riemannian [12] and non-Riemannian 
[13] acoustic metrics. The basic difference is that here we use a magnetic 
effective metric instead an acoustic one. Effective metrics in MHD maybe 
constructed, in near future by making use of the Navier-Stokes equation of 
viscous flows with nonlinear magnetic fields term. More recently we have 
shown that the non-Riemannian structure called acoustic torsion may also 
exist in this same case. The paper is organised as follows : In the section 2 we 
present the magnetic topology generalisation to gravitational torsion loops. 
In this section we also deal with the non-Riemannian loop magnetic metric 
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constrained by the MHD equations of dynamo theory. Section 3 shows that 
the more simple Riemannian geometry can take care of the MHD equations 
by making use of a time- independent diagonal metric. In section 4 discussions 
and conclusions are presented. 

2 Magnetic-like topology of torsion loops 

The Letelier teleparallel torsion loop metric is given by equation 

ds 2 = (dt + B.dxf - dx (1) 

when this metric is used on the Cartan's structure equation of the differential 
forms a vector relationship similar to the magnetic field derivative appears 
as 

J = VxB (2) 

— * — * 

where J is the gravitational analog to the magnetic field while B is the 
gravitational analog of the vector potential. It is clear that in analogy of 

— * — * 

magnetism we here also possess a vector gauge field where B — > B + Ve. 
From expression (2) we note that the torsion curve is really closed since the 
torsion vector field is divergence-free : V.J = 0. Of course now one can 
easily show that by defining the gravitational helicity by 

H g = J B.Jd 3 x (3) 

By analogy with the Berger and Field work [1] one has for the change in H 
according to the gauge freedom transformation that 

AH g = J Ve.Jd 3 x = J V.(eJ)d 3 x = J eJ.ndS (4) 

Where the integration of the first two integrals is on total volume of space 
V otherwise the field lines would close outside V. The Stokes theorem has 
also been used to obtain the last integral. Note, however, that only when 
the torsion vector is orthogonal to the vector n the helicity is conserved. To 
try to remedy this situation in the next section we drop the teleparallelism 
condition i? M Q| g 7 = (here a, (3 represents the four-dimensional spacetime co- 
ordinates) and consider the computation of other torsion components where 
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the helicity of the metric appears explicitly. Although the spacetime metric 
(1) is invariant with respect to a general coordinate transformation as should 
be the line elements in Riemannian geometry, it is not invariant with respect 
to above gauge transformations, and the torsion loop metric becomes 

#00 = 1 (5) 

g 0i = Bi + die (6) 

9ij = Sij + BiBj (7) 

where here v, /i — 0, 1, 2, 3 and latin indices takes values from one to three. 
It is clear that this metric exhibit explicitely the gauge freedom scalar. From 
the metric components is easy to compute the following components of the 
Cartan torsion tensor 

T i0j = \[di9oj - djgoi] = ^[diBj - djBi] (8) 

This expression can be recast in a more ellegant form by writing it in vector 
form as 

^o, = \[VxB} 1 (9) 

and 

T m = (10) 

where the Einstein summation convention was used in this last expression. 
Finally the last component of Cartan torsion is 

= i[a j£( vx B -)<] (id 

where e kh is the Levi-Civita symbol. From this last expresion we not that 
by contracting the indices I — i we note that a new generalised definition of 
gravitational helicity can be obtained since 

J e tkj T kij d 3 x = ^J Ve.VxSd 3 x = H g (12) 

which shows that this new definition coincides with the old with the advan- 
tage that now the full torsion tensor is consider and not only the torsion 
vector part. Since the component T i0 j can also be expressed in terms of the 
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vector J one may express the gravitational helicity by yet another integral 



as 



/ 



e lij BiT i0j d 3 x = l - J [B.WxB]d 3 x = H g (13) 



In the next section we propose a solution for the problem of helicity in 
Riemann-Cartan spacetime. In this section we show that it is possible to 
show that gravitational helicity is conserved as long as we extend the space- 
time to a more general Riemann-Cartan one instead of the teleparallel space- 
time. It is easy to show that by considering the torsion loop metric (1) in 
differential forms notation 

ds 2 = (uj ) 2 - (u 1 ) 2 - (u 2 ) 2 - (to 3 ) 2 (14) 

where the basis one-forms are 

u° = (dt + B.dx) (15) 

and uj 1 = dx, u 2 = dy and u 3 = dz. Now a small perturbation of connection 
one-forms in the teleparallel case according to the formula 

w\ = -\e in lp [J ] m u; n (16) 

where the new torsion vector field now in RC spacetime Jo is given in terms 
of the old teleparallel vector field by 

j = j_VxS (17) 

Note from this expression that the vector torsion field Jo represents also a 
loop since 

V.J = V.J = (18) 

where we have used the fact that V.[VxB] = 0. Note also that now the new 
definition of helicity similar to the previous one is 

H RC g = J B.J d 3 x (19) 

From this definition of helicity , where J was simply replaced by Jo, one is 
able to obtain the new expression for the variation of the helicity by 

AH RC g = J Ve.J d 3 x = J V.(eJ )d 3 x = J eJ .ndS (20) 



Note that from expressions (20) and (19) it is possible to obtain 



— * — * — * 

n. J = n.J - n.[VxB] = 



(21) 



This expression and (20) together lead to the conservation of the helicity in 
Riemann-Cartan spacetime given by 



One must notice that the condition that leads to this conservation does not 
imply that the torsion vector J is now orthogonal to the static torsion loop 
plane, which solves the contradiction this would imply in teleparallel space- 
time. In the spirit of this sectiona non-Riemannian loop magnetic metric 
can be obtained from the Letelier metric (1), the only difference, however, is 
that now the vector B in the metric coeficient is a true magnetic field. This 
choice is very convenient since by makind use of Cartan's calculus of differen- 
tial forms Letelier relation (2) is identical to the magnetic equation where the 

— * 

torsion vector field J is now equivalent to a electric current. Therefore in this 
case our system is equivalent to a circular current carrying loop generating 
magnetic fields. In the next section makes use of another magnetic metric 
where now the metric is Riemannian and the MHD equations is consider to 
constrain the Cartan's equations. The advantage of considering B as a real 
magnetic field is that relation (2) then implies a dynamo generating magnetic 
fileld from Cartan torsion. This idea has been fully sustained by De Sabbata 
and Gasperini [14]and more recently by Opher and Wichoski [15]. 

3 Riemannian magnetic metrics in MHD 

In this section we shall consider the diagonal magnetic metric given by the 
following line element in cylindrical coordinates adequate to treat the geom- 
etry of tubes in plasma physics 




(22) 



rfs 2 = ^2 _ B 2 dr 2 _ Bg 2 r 2 d9 2 _ ^2^2 



(23) 



where the basis one-forms are given by 



L0° = dt 



(24) 
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uj 1 = B r dr (25) 

u? = B e rde (26) 

uj 3 = B z dz (27) 
By making use of the first Cartan's structure equation 

T a = dco a + u a f3AojP (28) 

where the symbol A means the exterior product and u>^ 1 represents the con- 
nection one-form and T 13 represents the Cartan torsion two-form. All com- 
ponents of the magnetic field depends only upon the radial coordinate r. 
Since the metric is Riemannian the torsion forms vanish and this can be 
used together with the MHD equations for helical fields to constrain the ge- 
ometry and to find out the Riemann curvature of the magnetic metrics. To 
able to accomplish this task we consider the following MHD equations in the 
steady-state case leading to the phenomelogical Maxwell equation 

kvxB = VxB (29) 

where k is a constant. By considering the helical flow [11] 

B — (B r , Bg, B z ) (30) 

v = (0,ujr,v) (31) 

where u and v are constants. These vectors used in equation (29) yields the 
following conditions 

urB z = vB e (32) 
-Bj = kvB r (33) 
rBe = -kur 2 B r (34) 
Here the upper prime represents derivation with respect to the radial coor- 

— * 

dinate r. By making use of the solenoid condition V.B = one obtains the 
following solution of the helical MHD flow 

B r = — (35) 
r 

B e = - -kujcxr (36) 



7 



B z = C3 — kvc\lnr (37) 

To apply this MHD solution to constrain our Riemannian magnetic geometry 
one needs before to substitute the basis one forms into the Cartan equation 
(28) which yields 

T* = J r Auj r + ujtgAu + J z Auj z (38) 
T e = dJ + J r Au r + u e z Au z + uj\nJ (39) 
T z = u z r Au r + uJ z e Au e + uj z t AJ (40) 
T r = u r z Au z + u r e Au e + uj r t Auj l (41) 
where the only nonvanishing exterior derivative of the basis one-form is 

dJ = [rB'o + Be\drAd6 (42) 

Substitution of (42) into (39) yields the only nonvanishing component of the 
connection one-form 

J r = ^ ± ftl* (43) 

where q with (i = 1,2,3) are the integration constants. Substitution of 
Maxwell equation (34) into expression (43) yields the following constraint of 
MHD equations to the Riemannian geometry of magnetic metric 

J r = -kur 2 d6 (44) 

which from the second Cartan's structure equation 

R a p = du a p + ^^Au'p (45) 

yields 

R\ = dJ r (46) 
which along with the definition of the curvature two-form 

R a p = R a PlS ^Auj 5 (47) 

yields the following component for the Riemann tensor of the magnetic man- 
ifold 

E?rO Z = -2koor (48) 
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which is equivalent to the expression for the Riemann curvature of the mag- 
netic manifold 

R e rdz = -2kv(r) (49) 

This is the Riemannian curvature of the helical MHD flow. A similar relation 
between the curvature of intratube and the velocity of the flow has been 
obtained previously by Pelz [16] in the context of vortex filament models. A 
particular case of the above Riemannian MHD metric maybe consider as 

ds 2 = dt 2 - dx 2 - B y 2 dy 2 - dz 2 (50) 

where we have consider now the Cartesian retangular coordinates (x, y, z) 
for the spatial part of the magnetic metric. Here the only nonvanishing 
component of the magnetic field is given by B y (x,y). By again making use 
of Cartan's calculus of differential forms yields the following equations 

u 3 2 = [VxB] z dx (51) 

which in turn yields the following Bianchi identity 

dR 3 2 = duo 3 2 = d y d z [V xB] z dyAdxAdz (52) 

By considering equation (29) one obtains the equation 

dR 3 2 = du 3 2 = d y d z [kvxB] z dyAdxAdz = (53) 

which is consistent with the Bianchi identity dR a p = where d is the exterior 
derivative. In all the above example we found out an interesting interplay 
between the equations of plasma physics and the geometrical equations of 
the Cartan's calculus of differential forms. 

4 conclusions 

A natural extension of the magnetic topology of torsion loops to investigate 
knots in MHD can be undertaken by generalising the static Letelier's tor- 
sion loops to time-dependent torsion loops. This extension unfortunately 
has been proved very difficult [17] even by modern computation techniques 
using the OrtoCartan program. Non-Riemannian geometry and topology 
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of torsion curves is discussed. A new definition of the magnetic- like grav- 
itational helicity is proposed. We show that the extension of teleparallel 
spacetime to Riemann-Cartan spacetime allows us to possible applications 
of the mathematics discussed here in astrophysical models may be proposed 
in near future. Solitonic equations in non-Riemannian background may also 
be considered in near future. 
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